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ABSTRACT 

High resolution spectra from the Hinode EUV Imaging Spectrometer (EIS) have revealed that 
coronal spectral line profiles are sometimes asymmetric, with a faint enhancement in the blue wing 
on the order of 100 km s _1 . These asymmetries could be important since they may be subtle, yet 
diagnostically useful signatures of coronal heating or solar wind acceleration processes. It has also been 
suggested that they are signatures of chromospheric jets supplying mass and energy to the corona. 
Until now, however, there have been no studies of the physical properties of the plasma producing 
the asymmetries. Here we identify regions of asymmetric profiles in the outflows of AR 10978 using 
an asymmetric Gaussian function and extract the intensities of the faint component using multiple 
Gaussian fits. We then derive the temperature structure and chemical composition of the plasma 
producing the asymmetries. We find that the asymmetries are dependent on temperature, and are 
clearer and stronger in coronal lines. The temperature distribution peaks around 1.4-1.8 MK with 
an emission measure at least an order of magnitude larger than that at 0.6MK. The first ionization 
potential bias is found to be 3-5, implying that the high speed component of the outflows may also 
contribute to the slow speed wind. Observations and models indicate that it takes time for plasma 
to evolve to a coronal composition, suggesting that the material is trapped on closed loops before 
escaping, perhaps by interchange reconnection. The results, therefore, identify the plasma producing 
the asymmetries as having a coronal origin. 
Subject headings: Sun: corona — Sun: abundances — solar wind 



1. INTRODUCTION 

Understanding how magnetic energy is generated, 
transported, and dissipated to heat the solar corona, 
where the solar wind originates and how it is accelerated, 
remain the most important unsolve d problems in solar 
astrophysics. The Hinode mission (jKosugi et all 12007ft 
was designed to investigate these topics and its instru- 
ments are providin g many insights. Th e EUV Imaging 
Spectrometer fEIS. iCulhane et al~ll2007j ) is able to mea- 
sure the solar spectrum in the 171-212 A and 245-291 A 
wavelength ranges and provides measurements that can 
be used to derive plasma properties such as density, tem- 
perature, emission measure (EM), and chemical compo- 
sition. An important advance made possible by EIS is 
the ability to measure coronal line profiles at high spec- 
tral resolution (22 mA) . These measurements have shown 
that t he line profiles of ten have faint blue wing asymme- 
tries (jHaraieLaDIMQl). 

The origin and importance of these asymmetries has 
not been clearly established. Nanoflare reconnection 
models predict that ver y hot lines should show weak 
blue wing enhan cements (jPatsourakos fc Klimchukll2006t 
iKlimchukl l2006f l . while competing chromospheric jet 
models sugest that upflows sh ould be seen at loop foot - 
points at 'warm' temperatures (fAschwan den et al.ll2007D . 
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Active region outflows (ISakao et aLI l2007t lHarra et al.l 
[2008t iDel Zannal[2008tlDoschek et al.ll2008t) clea rly show 
asymmetries ( Brvans et al.ll2010t lDoschekll2012t) . so they 
could have some connection to the solar wind. They 
are less easy to associate w ith specific fea tures in other 
areas of an active region (jDoschekl I2012I ). but several 
studies have indicated that they may be s een in and 
around plage regions and loop footpoints (lHara et aLI 
120081: iDe Pontieu et al .l l2009l 120111; lPeterll2010D. thoug h 
they may be transient (|Ugarte-Urra fe Warrenl 1201 lh . 
Ther e are many ongoing investigations using EIS 
data dMcIntosh fc De Pontieull2009allbl: iTian et al.ll2011b 
iNishizuka fc Haral 120111: iMartmez-Svkora et al.l l201ltT . 
One difficulty is that the asymmetries, even when 
present, account for less than ~ 20% of the total area of a 
line profile, so multiple Gaussian fitting or novel analysis 
techniques have been developed to ex tract as much in- 
formation as possible from the profile s (jDe Pontieu et ail 
l2lMlBn^s et alJf20Tot |Pet^f20Toh . 

To date, studies of the asymmetries have focused on 
their location and velocity characteristics, but there have 
been no studies of the physical properties (temperature, 
chemical composition) of the plasma producing them. 
These diagnostics could provide new clues for coronal 
heating, solar wind, or chromosphere-corona connection 
studies, and this is the purpose of the investigation pre- 
sented here. We focus on the outflow regions of AR 
10978 that was observed by EIS in December 2007. AR 
10978 was chosen to assess the importance of the asym- 
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Fig. 1. — EIS Fe XII 192.394 A image of AR 10978 (left) and Doppler velocity map (right). The area used for analysis of the outflow is 
shown by the large white box. The size of the box is 28" X 21". Two smaller boxes show other areas discussed later. Contours indicating 
the locations where there is a > 10% blue-wing asymmetry in the line profile are also shown. Contours shorter than 2% of the largest one 
are omitted for clarity. The intensity image was treated with a Gaussian smoothing filter for presentation. 



metric component since the outflows of this region have 
previously been linke d to the solar wind thro ugh abun- 
dance measurements (|Brooks fe Warren 2011'). Further- 
more, strong blue wing asymmetries are clearly seen 
(| Br vans et alj|2010f ), enabling a more robust analysis in 
what is a difficult measurement. We use an asymmetric 
Gaussian function to detect deviations from a symmetric 
profile and then extract the intensities of the asymmet- 
ric component for a range of lines from Si X, S X, and 
Fe VIII-Fe XVII using multiple component fitting. We 
use these measurements to derive temperatures, emission 
measures, and relative abundances for the asymmetric 
component. 

2. DATA REDUCTION AND ANALYSIS 

We analyze observations of AR 10978 obtained on De- 
cember 10-15, 2007. The data were processed using stan- 
dard software. Context intensity and Dopppler velocity 
images are shown in Figure [TJ A large area of outflow is 
seen on the Western side of the AR. 

To detect regions of non-Gaussian shaped line profiles, 
we fit the data with an asymmetric Gaussian function of 
the form 



J(A) - 6(A) 




2a 2 1 A - An 

(A — An) 2 
2<rl 



A > An 



(1) 



where /(A) is the intensity at wavelength position A, An 
is the centroid of the asymmetric Gaussian, A is the area 
under the curve, ax and an are the widths of the Gaus- 
sian functions that reproduce the left and right wings 
of the profile, and 6(A) is a polynomial representing the 
background which is assumed to be linear. The asymme- 
try a is defined as (<jl — <tr)/((Tl + &r)- This function 
collapses to a Gaussian if there is no asymmetry Con- 
tours where a > 0.1 (greater than 10% asymmetry) are 
overlaid on the intensity and velocity maps in Figure [TJ 
They clearly show asymmetries in the outflows. Only a 
small fraction of the pixels in this raster scan (~ 9%) 
show line profile asymmetries as large as this, and only 
~ 2% show asymmetries greater than 20%. It is clear 



that the asymmetries are not very prevalent, at least in 
this region. 

To increase the signal for fitting the weaker lines, we 
average the spectra over boxed regions within the con- 
tours. An example box within the Western outflow on 
December 10 is shown in Figure [TJ To enable compari- 
son of the results with our previous work on this region, 
we first performed single Gaussian fits to the averaged 
spectra and extracted the intensities of a subset of Si X, 
S X, and Fe VIII-Fe XVII lines. 

To extract the intensities of the asymm etric compo - 
nents we performed multiple Gaussian fits. iPeterl (|2010l ) 
pointed out that fitting multiple Gaussians implicitly as- 
sumes that the emission is the combined contribution of 
a fixed number of spatial components. Thi s may not be 
true o f the real emission, as also noted by iBrvans et al.1 
(2010), which could be produced by a continuous dis- 
tribution of plasma velocities. It is important, there- 
fore, to consider uncertainties in the fitting procedure 
we adopt. We have found that fitting two components 
of equal wi dth leads to consis tent and reasonable results 
in general. iTian et all (|2011l ) have also found that the 
minor component usually has a Gaussian width compa- 
rable to the main component. There are circumstances, 
however, where this assumption is unsatisfactory: if the 
asymmetric component has a weak high speed tail, the 
restriction of fixed widths can lead to an underestimation 
of the minor component intensity. This sometimes affects 
the critical Si X 258.375 A line but not S X 264.233 A. 
Since S X 264.233 A is weaker, the high speed tail may 
be lost in the noise on these occasions, and we prefer to 
make a direct comparison by applying the same method 
to both lines. When S X 264.233 A is strong enough, 
however, we assume that any high speed tail should be 
visible, and we allow the width of the minor component 
to vary from the major component. In summary, in most 
cases we fit two components of fixed width to all the lines. 
In a few cases (3), however, we allow the minor compo- 
nent width of the Si X 258.375 A and S X 264.233 A lines 
to vary. Even in these cases all the Fe lines are fitted 
with fixed widths. Examples for spectral lines covering 
a range of temperatures are shown in Figure [U 
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Fig. 2. — Multiple component fits to the spectral line profiles for the boxed region of the Eastern outflow shown in Figure Q] Fits to 
Fe VIII 185.213 A, Si X 258.375 A, S X 264.233 A, Fe XII 192.394 A, Fe XIII 202.044 A, and Fe XV 284.160 A are shown, spanning a range of 
temperatures from 0.4-2.2MK. The histograms and green stars show the observed spectra and the solid green lines show the composite fit 
to the data. The red line shows the primary component of the fit, and the blue line shows the minor component (blue wing asymmetry). 



Once the intensities are extracted we perform the 
EM calculations. Recent analysis has shown that the 
outflows are more com plex than previous ly thought. 
I Warren et all (|2011[ ) and lYoung et al.l (|2012D have found 
that lower temperature lines show downfiows on loop 
structures that are unassociated with the upflows. To 
account for this, we only use lines where the total out- 
flow emission is blue-shifted. If a lower temperature line 
is red-shifted, the intensity is set to zero and a 3<r error 
is used for the uncertainty, where a is the error in the 
measured background intensity from the fit. This cor- 
rection is also applied to any line that is too weak to 
obtain a reliable fit. All the other lines have t he instru- 
ment calibration uncertainty (|Lang et al.ll2006f ) added in 
quadrature to the intensity measurement error. 

To derive the EM we use the Markov-Chain Monte 
Carlo (MCMC) algorithm dist ributed with the PINTo- 
fALE spectroscopy package (|Kashyap fe Drake! I1998L 
2000) and compute cont ribution functions from the 
CHIA NTI v7.0 database (|Dere et all 119971 ; lLandi et al.l 
I2012fh We adopt the CHIANTI ioniz ation fractions 
and t he photospheric abundances of iGrevesse et al.l 
(|2007f h The contribution functions are calculated at 
the electron density derived using the reliable Fe XIII 
202.044A/203.826A ratio for the total emission. The 
Fe XIII 203.826 A line is blended on the blue side, how- 
ever, with Fe XII 203.734 A, making it difficult to ex- 



tract any asymmetric component with reasonable uncer- 
tainties. Therefore, we use the same density to calcu- 
late the contribution function for the asymmetric com- 
ponent. We explored the possibility of using the Fe XII 
186.880 A/192. 394 A ratio instead and were able to ex- 
tract the asymmetric component for the 186.880 A line, 
but the known discrepa ncy between densiti es derived us- 
ing Fe XII and Fe XIII (|Young et al.ll2009f ) also became 
apparent for these outflows and introduced uncertainties 
around the temperature of formation of the critical S X 
and Si X lines used for the abundance measurements. 
For direct comparison between the asymmetric compo- 
nent and the bulk outflow, and also for benchma rking 
against our previous work (jBrooks fe Warrenl 120111 ). we 
decided to only use Fe XIII. 

The methodology used to derive the first ioniza- 
tion potential (FTP) bias has been presented by 
iBrooks fe Warren] (|201lD . and we have found our cho- 
sen lines to be reliable in separate emission measure 
analy s es of a variety of sol ar regi o ns dWarren fe Brooks! 
\2XM . IBrooks et all \2(M \2mi I Warren et al.l 1201 11) . 
Many of the m have also been independently verified in 
other studies <|Testa et al.ll2011t IWinebarger et al.ll2011l : 
iTripathi et alj|2011l ), strengthe ning our confidence in the 
atomic data. Furthermore, in IBrooks fe Warrenl ([201 If ) 
we tested the method by applying the analysis to eight 
polar coronal hole observations and verified that photo- 
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Fig. 3. — Emission measure results for all the emission in the outflow (left) and for the asymmetric component only (right). The blue 
line shows the best fit solution, la error bars computed from 100 Monte Carlo simulations using perturbed intensities are indicated. These 
show the possible spread in results. The green lines are emission measure loci curves indicating where the Fe lines constrain the solution. 
The dashed curves denote lines where the computed intensity was used (because the observed intensity was set to zero; see text). 



TABLE 1 

Differential Emission Measure Model 





Total Outflow Emission 


Asymmetric Component 




Line ID 


Tjfjs 0\T Idem R(%) 


lofts °\T Idem R(%) 


P(%) 



Fe VIII 185.213 


0.0 


5.3 


10.8 


0.0 


0.0 


5.9 


1.2 


0.0 


0.0 


Fe IX 197.862 


32.1 


7.1 


21.1 


34.3 


2.4 


0.6 


2.6 


6.9 


7.4 


Fe X 184.536 


161.4 


35.5 


131.6 


18.5 


22.6 


5.6 


18.2 


19.7 


14.0 


Fe XI 188.216 


278.7 


61.3 


313.4 


12.5 


42.2 


9.3 


45.9 


8.9 


15.1 


Fe XII 192.394 


150.4 


33.1 


177.7 


18.1 


23.8 


5.3 


24.4 


2.2 


15.8 


Fe XIII 202.044 


496.6 


109.3 


378.7 


23.8 


56.0 


12.4 


46.8 


16.5 


11.3 


Fe XIII 203.826 


195.1 


43.0 


148.7 


23.8 












Fe XIV 264.787 


99.0 


21.8 


107.0 


8.1 


12.5 


2.8 


12.6 


0.7 


12.6 


Fe XV 284.160 


609.4 


134.1 


655.1 


7.5 


73.1 


16.1 


76.7 


5.0 


12.0 


Fe XVI 262.984 


25.0 


5.5 


24.6 


1.7 


2.7 


0.8 


2.6 


3.6 


10.8 


Fe XVII 254.87 


0.0 


7.5 


1.2 


0.0 


0.0 


8.9 


0.1 


0.0 


0.0 


S X 264.233 


18.5 


0.6 






4.2 


17.3 






22.8 


Si X 258.375 


99.6 


1.1 






21.5 


0.9 






21.6 



lobs is the observed intensity, aj is the intensity error with the photometric calibration uncer- 
tainty added in quadrature. Idem is the calculated intensity. R is the absolute difference between 
the observed and calculated intensity. P is the ratio of the intensity of the asymmetric component 
to that of the total outflow emission. R and P arc expressed as percentages. 



spheric abundances were obtained, as expected for the 
presumed source of the fast solar wind. Briefly, the 
S X 264.233 A and Si X 258.375 lines are close in both 
wavelength and formation temperature, so the ratio is a 
potentially useful diagnostic for measuring abundances. 
Unfortunately, the ratio of these lines has some temper- 
ature and density sensitivity. To account for the tem- 
perature sensitivity we convolve the ratio with the EM 
distribution derived for the outflows. The EM is recon- 
structed from the Fe lines only. This minimizes uncer- 
tainties due to the choice of elemental abundances. The 
EM is then scaled to match the Si X 258.375 A intensity. 
This scaling is introduced to account for uncertainty in 
the Fe /Si abundance ratio and is generally small (less 
than 40% in the majority of cases). The intensity of the 
S X 264.233 A line is then predicted. Since S is a high 
FIP element and Si is a low FIP element, the ratio of 
predicted to observed S X 264.233 A intensity is the FIP 
bias. The inferred FIP bias is not very sensitive to the 
assumed density. Changing the density by a factor of 3 



changes the derived FIP bias by approximately 40%. 
3. RESULTS 

It is evident from Figure [5] that the magnitude of the 
asymmetric component is strongly dependent on tem- 
perature. Compare, for example, the fits to Fe VIII 
185.213 A and Fe XV 284.160 A. Table Q] shows both the 
total intensity and that of the asymmetric component 
for all the lines used for the EM analysis of the region in 
Figure [1] For this example, the asymmetric component 
amounts to 7% of the total emission at 0.7 MK, increasing 
to 16% at 1.6 MK. This result is broadly representative 
of all the regions we have studied. 

The derived EM distributions are shown in Figure |31 
The left panel shows the EM for the whole outflow at a 
derived density of log (n e /cm~ 3 ) = 8.55, and the right 
panel shows the EM for the plasma producing the asym- 
metries. Table Q] shows that the observed intensities are 
reproduced by the EM model to within 30% for all lines 
except Fe IX 197.862 A which is a little too bright in 
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TABLE 2 

Outflow Properties 





Total Outflow Emission 


Asymmetric Component 


Data 


Outflow T p EM P FIP 


T p EMp FIP AlogEM 



eis_ll_20071210-0019 


West 


1.6 


27.6 


3.3 


1.4 


26.7 


2.9 


3.4 


eisJl_2007121CL0019 


West 


1.8 


28.1 


3.5 


1.4 


26.9 


3.3 


2.5 


eisJl.20071210.0019 


West 


1.4 


27.6 


3.5 


1.8 


26.6 


3.4 


2.8 


eis_ll_20071211_1025 


West 


3.2 


27.8 


3.3 


1.6 


26.8 


3.0 


2.7 


eis_ll_20071211_1025 


East 


3.2 


28.3 


1.1 


1.6 


27.1 


4.7 


4.3 


eisJl.20071211.1025 


West 


1.8 


28.1 


3.2 


2.0 


27.0 


3.3 


2.8 


eisjl.20071211.1025 


West 


1.8 


27.5 


3.3 


1.4 


26.6 


2.8 


2.1 


eis_ll_20071211_1025 


West 


1.8 


27.8 


3.3 


1.6 


26.6 


3.8 


3.2 


eisJl.20071212.1143 


West 


1.6 


27.6 


2.9 


1.6 


26.7 


3.0 


3.2 


eis_ll_20071212_1143 


East 


1.6 


28.0 


4.6 


1.4 


26.9 


4.9 


3.7 


eis_ll_20071212_1143 


West 


1.8 


28.0 


3.6 


1.6 


26.8 


3.1 


2.4 


eis_ll_20071212_1143 


East 


2.2 


27.9 


3.4 


1.8 


26.7 


3.0 


4.1 


eisjl.20071213.1218 


West 


3.2 


28.2 


2.9 


1.6 


26.7 


4.9 


4.1 


eis_ll_20071213.1218 


East 


1.6 


27.7 


1.0 


1.4 


26.8 


3.8 


2.1 


eis_ll_20071215.0013 


East 


1.8 


27.8 


3.1 


1.4 


26.9 


4.2 


3.2 



The EIS datascts used arc given along with the location of the outflow. The peak tem- 
perature (in MK), peak emission measure (cm -3 ), and FIP bias are given for calculations 
using both the total emission and the asymmetric component only. The ratio of the peak 
emission measure to that measured at 0.6MK is also shown. 



the bulk emission case. This is seen in another example 
and in fact more often in Fe X 184.536 A. As discussed 
above, the transition from downflows on closed structures 
to outflows on open held lines takes place in this temper- 
ature range, and our method of abruptly including or 
excluding lines based on Doppler shift probably does not 
account perfectly for the gradual transition from emis- 
sion dominated by one type of structure to the next. In 
fact it is difficult to select boxes on the Eastern side of 
the AR that are both within the outflow and asymmetry 
regions and are clear of the fans. This is the reason that 
we have selected relatively more areas from the West- 
ern side in the following analysis of a sample of outflows. 
Note, however, that there is an Fe XI line at 184.412 A 
that is easily separable from Fe X 184.536 A in normal 
solar conditions, but becomes blended with the profile 
when there is a high speed component. This could affect 
the Fe X intensity measurements. Also, there is some un- 
certainty as to whether the Fe VIII and Fe IX intensities 
can simultaneously b e reproduced by the EM of every 
type of solar feature (|Schmeldl2012l ). 

To study the generality of the results, we repeated our 
analysis for 14 additional regions in these AR observa- 
tions. The chosen regions range in size from 10" x 8" to 
28" x 16". Our results, which are summarized in Tabled 
show that the temperature distributions for the outflows 
peak in the range 1.4-3.2 MK with an emission measure 
of log (EM/cm -5 ) = 27.5-28.3. The measured densities 
(not shown) are log (n e /cm -3 ) = 8.4-8.9. These num- 
bers a re consistent with previous EM analysis of AR out- 
flows dBrooks fe Warrenll201lHSlemzin et al.ll2012ft . The 
derived FIP bias falls in the range 2.9-4.6. 

The EM distributions for the asymmetric component 
appear to peak in a narrower temperature range: 1.4- 
2.0 MK, and they show rapidly decreasing emission mea- 
sure at temperatures below 1MK. We characterize this 
fall off by measuring the ratio of the EM at the peak to 
that at 0.6 MK. In all cases the EM is at least two or- 
ders of magnitude larger at the peak temperature. The 



derived FIP bias for the asymmetric component falls in 
the range 2.8-4.9, which is similar to that of the bulk 
outflow. 

Although we excluded lines with a bulk red-shifted pro- 
file, they often have weak blue wing emission, see e.g., 
Fe VIII 185.213 A in Fi gure[2] We therefore re-computed 
the results for the asymmetric component including these 
lines. We found that they are far enough away in tem- 
perature that the peak emission measure, temperature, 
and FIP bias are not affected significantly if they are 
included. In contrast, A log EM obviously changes sig- 
nificantly for individual cases, though, as discussed, it 
is unclear if the low temperature EM becomes too large 
when contaminated by emission from other structures. 
The range of results remains broadly the same, however, 
and in all cases the EM is still at least an order of mag- 
nitude larger at the peak temperature. 

4. SUMMARY AND DISCUSSION 

We have measured the physical properties of several 
outflow areas in AR 10978 during its disk passage in De- 
cember 2007. We also extracted measurements for the 
plasma producing the faint blue wing asymmetries in 
the EIS line profiles. Results for the bulk outflow are 
consistent with those of previous studies. 

The new measurements for the asymmetric compo- 
nent show that the magnitudes of the asymmetries are 
strongly dependent on temperature and peak in the 
corona. The emission of the source plasma appears also 
to be dominated by, and come from, a coronal tempera- 
ture range. The relative abundance measurements indi- 
cate a similar FIP fractionation value to that of the bulk 
outflow. 

These results have two important consequences. First, 
the FIP bias values (3-5) are slightly higher than but 
br oadly similar to those of the slow speed solar wind (3- 
4. Ivon Steiger et all 12000ft . The higher value probably 
reflects the difficulty of making this measurement, but 
indicates that the high speed component of the outflows 
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could also be a contributor and source of the slow speed 
wind. 

Second, in relation to coronal heating studies, pre- 
vious work has indicated that even if rapid fractiona- 
tion occurs locally, it takes at least several hours for the 
plasma to evolve to a coron al composition on a large scale 
(jFeldman fc Widine] I2003D . Other work has suggested 
that the material is stored on coronal loops and is later 
released by interchange reconnection with open field lines 
(I Wang et a.l.lll996L ISchwadron et al.lfl999h . Our FIP bias 
results therefore imply that the asymmetries in the AR 
outflow are likely to be produced by a mechanism that 
releases plasma of coronal origin. This is also consistent 
with our determination that the EM distributions are 
dominated by coronal emission. 

Finally, it is unclear whether our results could sup- 
port a picture whereby the asymmetries are a signa- 
ture of plasma that is rapidly injected into the corona 
from the chromosphere. In addition to the apparent 
coronal source, recent models suggest that fractionation 
is harder and therefore reduced for plas ma that passes 
rapid ly through the fractionation region (jLamin 
2012). In the most obvious case, we would expect to 



measure a photospheric composition, and therefore the 
S line would be relatively more asymmetric than the 
Si line; the opposite of what we see. We stress, how- 
ever, that this initial study has been of asymmetric line 
profiles averaged over relatively large areas in AR out- 
flows only. Similar studies on smaller scales may also 
be possible in bright moss or plage, for example, since 
good line profiles are obtained with EIS even with short 
exposures ( Brooks k. Warrer3l2009fl . Linking these weak 
asymmetries to small scale chromospheric jets may re- 
quire a higher sensitivity spectrometer with better spa- 
tial resolution covering the whole atmosphere from chro- 
mospher e to corona. Such a spectrometer is proposed for 
Solar-C (|Teriaca et al.l[20H . 

We thank Martin Laming for helpful discussions. This 
work was performed under contract with the Naval Re- 
search Laboratory and was funded by the NASA Hinode 
program. Hinode is a Japanese mission developed and 
launched by ISAS/JAXA, with NAOJ as domestic part- 
ner and NASA and STFC (UK) as international partners. 
It is operated by these agencies in co-operation with ESA 
and NSC (Norway). 
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